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Abstract. Through polarization observations, circumstellar masers are excellent probes of the magnetic field in the envelopes 
of late-type stars. Whereas observations of the polarization of the SiO masers close to the star and on the OH masers much 
further out were fairly commonplace, observations of the magnetic field strength in the intermediate density and temperature 
region where the 22 GHz H 2 masers occur have only recently become possible. Here we present the analysis of the circular 
polarization, due to Zeeman splitting, of the HtO masers around the Mira variable stars U Her and U Ori and the supergiant 
VX Sgr. We present an upper limit of the field around U Her that is lower but consistent with previous measurements, reflecting 
possible changes in the circumstellar envelope. The field strengths around U Ori and VX Sgr are shown to be of the order of 
several Gauss. Moreover, we show for the first time that large scale magnetic fields permeate the circumstellar envelopes of an 
evolved star; the polarization of the H2O masers around VX Sgr reveals a dipole field structure. We discuss the coupling of the 
magnetic field with the stellar outflow, as such fields could possibly be the cause of distinctly aspherical mass-loss. 
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1. Introduction 

At the end of their evolution, a large majority of stars go 
through a period of high mass loss while climbing the asymp- 
totic giant branch (AGB). This mass loss, of the order of 10~ 7 
to 10~ 4 M yr~', is responsible for the formation of a circum- 
stellar envelope (CSE). It is also the main source for replen- 
ishing interstellar space with processed material. AGB stars 
such as the Mira variable stars, with main sequence masses 
less than a few M , will eventually develop as planetary neb- 
ulae (PNe). The heavier evolved stars such as the supergiants 
will eventually explode as a supernova (SN). The exact role 
of magnetic fields in the mass loss mechanism and the for- 
mation of CSEs is still unclear but could be considerable. 
The study of several maser species found in CSEs has re- 
vealed important information about the strength and structure 
of magnetic fields throughout the envelopes surrounding the 
late-type stars. At distances from the central star of up to sev- 
eral thousands of AU, measurements of the Zeeman effect on 
OH masers indi cate ma g netic fields strengths of a few milli - 
Gauss (e.g. lSzvmczak & Cohen! 1 19971 iMasheder et allll999l) . 
Addition ally, weak alig nment with the CSE structure is found 
(e.g. Etoka & Diamond, 2004). Observations of SiO maser po- 
larization have shown highly ordered magnetic fields close 
to the central star, a t radii of 5-10 AU where the SiO maser 
emission occurs (e.g. lBarvainis. Mcintosh. & Predmorelll987t 
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iKemball & Diamon 3 Il997l) . When interpreting the circular 
polarization of the SiO masers as standard Zeeman splitting, 
the magnetic field strength determined from these observations 
could be up to several tens of Gauss. However, a non-Zeeman 
interpretation can explain the observations by magnetic field 
streng ths of only several tens of milliGauss dWiebe & Watsorl 
1998). 

Recently, high circular polarization of circumstellar 
H2O maser s was observed for a small sample o f late- 
type stars dVlemmings. Diamond. & van Langeveldel I2002L 
hereafter V02). H2O masers occur at intermediate dis- 
tances in the CSE, in gas that is a factor of 10 - 

100 more dense than the gas in which OH m asers oc- 

I n 1 1 1 

cur. Vlemmings. Diamond. & van Langevelde (2001, hereafter 

V01) managed to detect the Zeeman splitting of the circumstel- 
lar H2O maser around the supergiant star S Per, even though for 
typical field strengths of a few hundred mG the H2O Zeeman 
splitting is extremely small, only m 10 3 times the typical half- 
power width of the H2O maser line (Av^ » 30 kHz). Additional 
detections were discussed in V02 and typical field strength of 
between 0.2 and 1 G were found, usin g the full, no n-LTE radia- 
tive transfer equations from Nedoluha & Watsorj Jl992l) . Such 
high magnetic field strengths at a few hundred AU from the 
central star give strong support for the standard Zeeman inter- 
pretation of the SiO maser polarization. 

The results from V02 indicate that the magnetic field 
strengths measured in the circumstellar maser regions are con- 
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sistent with either a r~ 2 field strength dependence on the in- 
creasing distance to the star, similar to a solar-type magnetic 
field or possibly an r~ 3 dependence as produced by a dipole 
magnetic field. This implies surface magnetic fields of hun- 
dreds to several thousand Gauss, indicating that the magnetic 
field pressure dominates the thermal and radiation pressure at 
the surface of the star. As a result, the magnetic field likely 
plays a very important role in driving the late-type star mass 
loss. Additionally, magnetic fields possibly play an important 
role in shaping the CSEs. This will also affect the formation 
of PNe, that are often observed to have distinctly non-spherical 
shapes. 

Here we present new results on the magnetic field strength 
in the envelopes of several late-type stars. We discuss the obser- 
vations and data calibration in §|5]and the analysis method in 
§|3] The new and previous results for our sources are presented 
in §E In §[5]the specific results on the supergiant VX Sgr are 
discussed, where the maser polarization observations allow for 
the investigation of the morphology of the magnetic field. The 
consequences of these measurements are discussed in §|6]and 
are followed by our conclusions in §0 

2. Observations 

The observations were performed at the NRAO Very Long 
Baseline Array (VLBA) on April 20 2003. The average beam 
width is 0.5 x 0.5 mas at the frequency of the 6k, - 523 ro- 
tational transition of H2O, 22.235 GHz. We used 4 baseband 
filters of 1 MHz width, which were overlapped to get a veloc- 
ity coverage of » 44 km/s, covering most of the velocity range 
of the H2O masers. Similar to the observations in V02, the data 
were correlated multiple times. The initial correlation was per- 
formed with modest (7.8 kHz= 0.1 km s ) spectral resolution, 
which enabled us to generate all 4 polarization combinations 
(RR, LL, RL and LR). Two additional correlator runs were per- 
formed with high spectral resolution (1 .95 kHz= 0.027 km s ) 
which therefore only contained the two polarization combi- 
nations RR and LL. This was necessary to be able to detect 
the signature of the H2O Zeeman splitting in the circular po- 
larization data and to cover the entire velocity range of the 
H2O masers. Each source-calibrator pair was observed for 6 
hours. The calibrator was observed for 1.5 hours in a number 
of scans equally distributed over the 6 hours. 

2.1. Calibration 

The data analysis p ath is described in detail in V02. It fol- 
lows the method of iKemball. Diamond. & Cotton! Jl995l) and 
was performed in the Astronomical Image Processing Software 
package (AIPS). The calibration steps were performed on the 
data-set with modest spectral resolution. We were also forced 
to flag several channels that suffered from strong interference. 
Fringe fitting and self-calibration were performed on a strong 
maser feature. The calibration solutions were then copied and 
applied to the high spectral resolution data-set. Finally, cor- 
rections were made for instrumental feed polarization using 
a range of frequency channels on the maser source, in which 
the expected frequency averaged linear polarization is close to 



zero. Then image cubes were created for stokes I,Q,U and V in 
the modest spectral resolution data set, and for stokes I and V 
in the high spectral resolution data-set. 

2.2. Sources 

We observed 4 late-type stars, the supergiant VX Sgr and the 
Mira variable stars U Ori, R Cas and U Her. U Her was previ- 
ously observed in December 1998 as part of the observations 
described in V02. Unfortunately, we were unable to detect the 
H2O masers of R Cas, presumably due to source variability. In 
Table ^ we list the observed sources with type, position, dis- 
tance, period and velocity, as well as the peak H2O maser flux 
at the epoch of observation. In the high spectral resolution to- 
tal intensity channel maps, the noise is dominated by dynamic 
range effects and is w 60 mJy. In the circular polarization po- 
larization maps the rms noise is m 20 mJy. 

The polarization of the OH masers of our sources has 
previously been observed and SiO maser polarization mea- 
surements have been performed on VX Sgr. Polarization 
of the 1612 MHz OH masers a round U Ori ind icate a 
magnetic fie l d of £ OH ~ 10 mG dReid et all Il979h . while 
IPalen & Fix! fcOOOl) found a field of B 0H * 1 mG on the 
1665 and 1667 MHz OH masers around U Her. The mag- 
netic fields in the main-line and satellite-line OH maser 
regions around VX S gr we re measured using MERLIN 
dChapma n & Co hierl Il986t iT rigilio. Umana. & Cohenl 
Il998l ISzvmczak Cohen. & Richards! 1200 ll) and the VLA 
dZell & Fixlll996l) . The polarization of the 1665 and 1667 MHz 
OH masers indicated a field strength of Boh ~ 2 mG, while 
the field in the 1612 MHz maser region is Boh ~ 1 mG. 
In the single dish observations of the circular polarizatio n 
of SiO masers by Barvain is. Mcintosh. & Predmorel {l987), 
they observe a circular polarization percentage of 8.7%. This 
indicates B S io ~ 90 G for VX Sgr. 

3. Method 

For the analysis of the polarization spectra we used both the ba- 
sic LTE interpretation and the full radiative transfer non-LTE 
interpretation, which were thoroughly described in V02. The 
LTE analysis consists of a standard Zeeman interpretation as- 
suming LTE. As a result, the narrowing and rebroadening of the 
maser profile caused by saturation are not reproduced. We have 
included the possibility of multiple masing hyperfine lines. In 
this analysis, the magnetic field strengths are determined by 
fitting a synthetic V-spectrum, proportional to dl/dv, to the po- 
larization spectra. The magnetic field strengths follow from 

Py — (^max ^min)/^max 

= 2 • A F _F ■ fi [Gauss] cos6yAv L [km s" 1 ]. (1) 

with Py the percentage of circular polarization, V max and V m i n 
the maximum and minimum of the synthetic LTE V-spectrum 
fitted to the observations. 7 max and Avl are the peak flux and 
the full width half-maximum (FWHM) of the maser feature re- 
spectively. B is the magnetic field strength at an angle 8 from 
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Table 1. Star Sample 



Name 


Type 


RA (J2000) Dec (J2000) 


Distance 


Period 


v rad 


Peak flux 






(* ms ) C" ") 


(pc) 


(days) 


(km/s) 


(Jy) 


VX Sgr 


Supergiant 


18 08 04.0485 -22 13 26.614 


1700" 


732 


5.3 


43.2 


UHer 


Mira 


16 25 47.4713 +18 53 32.867 


m b 


406 


-14.5 


9.75 


UOri 


Mira 


05 55 49.1689 +20 10 30.687 


300 c 


368 


-38.1 


6.44 


RCas 


Mira 


23 58 24.8725 +51 23 19.703 


176* 


430 


26.0 




" Marvel 


fl996),''Vlemmin2s et al. (2003), r KnaDD et al. (1998) 



the line of sight. The Ap_p coefficient depends on the masing 
hyperfine component and in V02 was estimated to be a 0.01 1, 
which is used here for all LTE estimates. 

For the non-LTE analysis, the coupled equations of state 
for the 99 magnetic su bstates of the three d ominant hyperfine 
components from lNedoluha & Watsonl i 19921) were solved for 
a linear maser in the presence of a magnetic field. The mod- 
eled spectra were then directly fitted to the observed spectra. 
When a direct model fit was not possible, we used Eq^with 
Ap_p = 0.018, which is the corresponding best estimate for the 
coefficient in the non-LTE case. 

The spectral fitting for both the non-LTE as well as the LTE 
analysis requires the removal of the scaled down total power 
spectrum from the V-spectrum to correct for small residual gain 
errors between the right- and left-polarized antenna feeds. It 
was found in V01 and in V02 that generally, the V-spectra were 
narrower than expected in the LTE case. This could be partly 
explained by introducing the non-LTE method, however, some 
narrowing remained. In V02, we attributed this to the applica- 
tion of uni-directional maser propagation in the modeling, and 
we found that a bi-directional maser could explain some of the 
narrowing. However, for computational efficiency we continue 
to use the uni-directional method, while allowing for additional 
narrowing of the V-spectrum, as we found that the estimates for 
the fractional circular polarization are robust. 



4. Results 

We have determined the magnetic field strength on the 
strongest maser features in the band (44 km/s wide). The results 
of both the LTE and non-LTE analysis are shown in Table[2]for 
features with peak fluxes down to ~ 5% of the brightest maser 
spot. As seen in column 6, the percentage circular polarization 
is relatively large and varies between 0.3% and 20% with the 
highest polarization being detected for U Ori. Magnetic field 
strengths obtained using the LTE Zeeman method are given in 
column 7 while the results obtained using the non-LTE mod- 
els are shown in column 8. Features that did not show circular 
polarization at a level higher than 3 times the rms noise in the 
V-spectrum are considered non-detections. However, the for- 
mal errors on the magnetic field strength also include the errors 
on total intensity and line width, and are thus larger, resulting in 
occasional magnetic field strengths around the lcr noise level. 
For the non-detections we have determined the absolute upper 
limits using the Ap_p coefficient obtained in the non-LTE ra- 
diative transfer models. In the LTE analysis the limits would be 
larger by a factor of « 1 .4 



Maser features that had circular polarization spectra that 
could not be fitted accurately are labeled in Table. |2] This 
is likely due to blending of several features, both spatially 
and along the velocity axis. The values for the magnetic field 
strength for those features were determined using only the 
minimum V m ; n and maximum V max of the V-spectrum after 
it was forced to be point symmetric. The features labeled ** 
(VX Sgr n, shown in Fig@]and VX Sgr w) suffered from exter- 
nal interference, most of which we managed to remove. Some 
effects remain however, and it is unclear if this had an effect on 
the measured magnetic field strength for those features. 

As in V02, we were unable to detect any significant linear 
polarization above a limit of ~ 0.02% for the strongest maser 
features and a few % for most of the weakest. There was a hint 
of linear polarization on the strongest maser feature of U Ori, 
but as the level of polarization was at most only at the lcr level 
of as 2% it cannot unambigiously be classified as a detection. 

4.1. UHer 

After the observations of U Her in V02 revealed a magnetic 
field in the H2O maser region of ~ 1.5 G, U Her was ob- 
served again in our second observational run. The integrated 
total intensity map of the H2O maser features is shown in 
Fig. n Even though the strength of the U Her maser features 
is similar in both observational epochs, we did not detect any 
significant circular polarization in the observations presented 
here. The upper limit we find is ~ 130 mG on the strongest 
and ~ 600 mG on the weakest maser feature. However, the 
spectrum and spatial distribution of the masers has changed 
significantly in the more than 4 years between the observa- 
tional epoch of December 13th 1998 and April 20th 2003. 
The velocity of the maser features observed in V02 (-19.3 
- -17.6 km/s) was several km/s more blue-shifted than the 
maser features observed here, that are close to the stellar ra- 
dial velocity of -14.5 km/s. Whereas the spatial distribution of 
the maser spots in V02 showed an elongated structure in the 
N-S direction, Fig. ^ indicates a significantly different struc- 
ture. The weakest VLBI feature detected at -15.82 km/s (la- 
beled c in Fig. Q corresponds in velocity to the feature that 
was found to be the strongest in t he MERLIN observations of 
Iviemmings. van La ngevelde. & Diamond (2002). In those ob- 
servations, this feature was found to be aligned with the optical 
position to within 15 mas and interpreted as the amplified stel- 
lar image. Assuming that the current feature at -15.82 km/s is 
also the amplified stellar image we obtain a lower limit on the 
distance on the sky between the brighter H2O maser features 



4 

Table 2. Results 



W.H.T. Vlemmings et al.: Magnetic Fields around Late-type stars 



Name 


Feature 


Flux (I) 


Vrad 


Av L 


Pv 


B cose LTE 


B cosff non-LTE 






(Jy) 


(km/s) 


(km/s) 


(xlO- 3 ) 


(mG) 


(mG) 


VX Sgr 


a 


18.8 


-2.01 


0.47 


21.0 


350±72 


512±90 




b* 


11.7 


-1.28 


0.38 


22.6 


-307±99 


-248±62 




c 


6.9 


7.18 


0.61 


171.2 


3705 ±249 


3375±265 




d 


23.1 


7.16 


0.52 


32.5 


-609±67 


-469±60 




e 


18.9 


6.81 


0.53 


19.2 


359±79 


249±46 




f 


36.7 


16.43 


0.57 






< 54 




g 


43.2 


13.77 


0.58 


3.5 


-75±34 


-63+20 




h 


9.6 


11.86 


0.57 


15.5 


-326±184 


-290+115 




i 


23.3 


12.73 


0.48 






<99 




j 


9.4 


-0.9 


0.62 


16.4 


-370±202 


-325±126 




k 


6.3 


-0.7 


0.56 


14.4 


285+281 


277±195 




1 


14.9 


6.73 


0.51 


66.7 


1231+123 


966±120 




m 


9.7 


20.6 


0.52 


9.9 


-188±140 


-167±87 




n** 


14.5 


-0.1 


0.87 


162.7 


-4954±238 


-4082±305 




o 


6.8 


21.19 


0.71 






< 381 




P 


6.3 


22.17 


0.62 






<471 




q 


8.5 


20.97 


0.49 






<239 




r 


3.4 


20.93 


0.62 






<687 




s 


2.7 


9.95 


0.54 






<987 




t 


8.4 


16.26 


0.48 






< 217 




u 


9.9 


13.64 


0.58 


4.5 


-94+154 


-85+99 




V 


7.2 


7.1 


0.51 


39.9 


-738±194 


-366+75 




w" 


4.7 


1.17 


0.56 


109.6 


2160+387 


2195+251 




X 


5.1 


18.24 


0.69 






< 512 




y 


6.8 


2.52 


0.46 


61.1 


1012+281 


921 ±208 


UHer 


a 


9.75 


-14.56 


0.55 






< 127 




b 


3.75 


-14.42 


0.58 






< 418 




c 


2.09 


-15.82 


0.55 






<610 


UOri 


a* 


6.44 


-41.05 


0.61 


22.6 


-517±152 


-401 ±96 




b 


1.87 


-41.23 


0.75 


217.9 


-5770±718 


-3192±425 




c 


0.95 


-40.85 


0.66 


140.4 


3462±1210 


1398+382 



* no accurate fit possible (see text) 
** suffer from interference (see text) 



and the star. For a distance to U Her of 277 pc, the projected 
separation between the star and the two brightest maser fea- 
tures is a: 16 AU. Since the velocity of these maser features is 
close to the stellar velocity, we expect the true separation to be 
similar. Assuming that the maximum magnetic field strength in 
V02 was measured at the inner edg e of the H?Q mase r region, 
which was found to be « 12.5 AU (iBains et all 120031) . the ex- 
pected magnetic field strength at 16 AU is « 700 mG assuming 
a dipole magnetic field configuration. The difference between 
the expected field strength and the measured upper limits could 
be due to a large angle 6 > 75° between the magnetic field 
lines and the line of sight. Alternatively, some of the magnetic 
field strength difference could be due to different maser clump 
density by up to a factor of ~ 10, if the magnetic field lines 
are frozen into the dense clumps. We could also be observing 
actual variability in the magnetic field. 

4.2. U Ori 

Fig- HI shows the integrated total intensity H2O maser map of 
U Ori. Very little structure is observed and the maser features 
seem to be somewhat blended both spatially as well as in ve- 



locity. Strong circular polarization up to =s 22% is detected on 
the three bright central features where interestingly the fea- 
tures b and c, to the left and right of the brightest feature a, 
show an opposite sign for the magnetic field. The total inten- 
sity and circular polarization spectra for these features are also 
shown in Fig.|2] From these we estimate the absolute magnetic 
field strength in the maser region to be ~ 3.5 G. While the 
central feature a does not show such a high field strength, the 
feature seems to be a blend of two features with oppositely di- 
rected magnetic fields. The more complex circular polarization 
spectrum of U Ori a can be explained by two features with a 
magnetic field of ~ 800 mG of opposite sign. The high mag- 
netic field strength measured for the Mira variable star U Ori is 
consisten t with the high field (10 mG) measured using the OH 
masers jReid et.allll97^> assuming a dipole r 3 dependence of 
the magnetic field on the distance to the star. This assumes that 
the H2O masers occur in a shell wit h an outer edge of 30 AU, 
as determined bv lBains et ai1(l2003l) . 
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MilUARC SEC i n 1 a in 

Cl-iii^i' :u RA !)5 49. KjNWIXI DLC 20 10 <I).<,K7000 1 U 

Peak flux = 6.3082E+0OJY/BEAM 
Levs=6.308E-01* (0.500, 1,2,4,6,9) 

Av [km/sj Av [km/sj 



Fig. 2. (left) Total intensity image of the H2O maser features around U Ori. (right) Total power (I) and V-spectra for selected 
maser features of U Ori. The bottom panel shows the best fitting synthetic V-spectrum produced by the standard LTE Zeeman 
interpretation (dashed line). The middle panel shows the best non-LTE model fit (dotted line). The corresponding total power fits 
are shown in the top panel. The V-spectra in the lower two panels are adjusted by removing a scaled down version of the total 
power spectrum as indicated in §|5] which is different for the LTE and non-LTE fits. The solid vertical lines show the expected 
position of the minimum and maximum of the V-spectrum in the general LTE interpretation. 



4.3. VXSgr 

The integrated total intensity H2O maser map of the supergiant 
VX Sgr is shown in Fig. [5] VX Sgr shows a rich structure of 
strong H2O maser features in a 350 x 250 mas ellipse, corre- 
sponding to 600 x 425 AU at a distance of 1 .7 kpc. The maser 
feature distri bution is similar t o that seen in the MERLIN ob- 
servation bv iMurakawa et all J2003L hereafter M03). Circular 
polarization was detected ranging from 0.3-16% for a large 
number of maser features as seen in Table [5] Several of the 
spectra are shown in Figs. [3] and |4] As mentioned above, a 
few of the maser features suffered from blending or external 
interference. The maximum magnetic field strength measured 
is ~ 4 G. Across the maser features we find a clear transition 
between a negative magnetic field in the S-E to a positive mag- 
netic field in the N-W. This is the first detection of a large scale 
magnetic field in a circumstellar H2O maser region. 

5. The Morphology of the Magnetic Field around 
VX Sgr 



Sgr has been studied exten- 
L 

1986| l7-ell & Fix! 119961 



The maser shell of VX 

sively with MERLIN and V LA ob se rvations (e .g. |Lane, 
1 19841 Chapman & Cob 



iTrigilio. Umana. & Cohen! Il998l M03). The recent ob- 
servations of M03 indicate that the H2O masers arise in a 
thick shell between 100 AU and 325 AU from the star. In this 
shell, the H2O masers seem to be accelerated from ~ 10 km/s 
to 20 km/s. The H2O maser shell shows a clear elliptical 
asymmetry which was modeled in M03 as a spheroidal maser 
distribution intersected with an under-dense bi-conical region 
at an inclination angle i = 60° + 30° from the li ne of sight and 
rotated = 200° + 30° on the plane of the skv. lMarvef (ll996) 
also observed the elliptical H2O maser distribution with the 
VLBA and, constrained by proper motions, fitted the emission 
region using an oblate spheroid with a maximum radius of 
250 mas at an inclination angle i ~ 46° and projected position 
angle ~ 230°. 

Observations of the 1612 MHz OH masers at ~ 1400 AU 
from the star also show the elliptical asymmetry observed 
in the H2O masers. The magnetic field strength in the 
1612 MHz OH maser r egion was estimated to be of th e 
order of a milliGauss (| Trigilio. Umana. & Cohenl ^998). 
ISzvmczak. Cohen. & Richards! J2001I) observed the linear po- 
larization of the OH masers and found structure that could be 
explained by a dipole magnetic field with an inclination axis 
of about 20° or 30°. Additionally, the observations of Zeeman 
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so 200 150 ioo 50 o -50 Av [km/s] Av [km/s] 

MilliARC SEC 
Center at RA 18 08 04.0485000 DEC -22 13 26.614000 
Peak flux = 4.61 12E+01 JY/BEAM 
Levs = 4.61 1E+00 * (0.500. 1, 3, 6, 9) 



Fig. 3. Similar to Fig.|2]for VX Sgr. 




-1 1 -1 1 -1.5 -1 -0.5 0.5 1 -1 -0.5 0.5 1 1.5 



Av [km/s] Av [km/s] Av [km/s] Av [km/s] 

Fig. 4. Similar to the spectra in Fig. |2] and Fig.[3]for additional maser features around VX Sgr. Feature VX Sgr n suffers from 
some interference effect in the right wing of the spectrum but this does not affect the magnetic field determination. 



pairs in the OH emission indicates a change in magnetic field 
direction across the maser region. This w as suggested to b e 
the signature of a dipole magnetic field JZell & FixL [l996). 
Further observations indicated that the magneti c axis was pro- 
jected at a position angle = 210°+ 20° JSzvmczak & Coherl 
I1997I) . Observations of the OH mainline masers emission at 
1667 MHz closer to the star are c onsistent with the same mag- 
netic angle (Richar ds et al! l2000). 



We have fitted a dipole magnetic field to our H2O maser 
magnetic field observations of VX Sgr. As the positions of 
the maser features in the shell along the line of sight are un- 
kno wn, we used th e acce lerating spherical outflow model of 
Chapman & Cohei] dl986l) to describe the velocity structure of 
the maser shell. Thus, the observed radial velocity of the maser 
features directly maps into the third spatial dimension, allow- 
ing for a three dimensional fit. A maximum likelihood fit was 
made to the dipole field, solving for the stellar position, the in- 
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Fig. 5. The best fitted dipole magnetic field for the H2O maser observations around VX Sgr (denoted by the solid circle), (left) 
The distribution of the H2O maser features indicating the measured magnetic field strengths. Open symbols denote a positive 
magnetic field while the closed symbols correspond to a negative magnetic field. The crosses represent the upper limits. The 
symbols have been scaled relative to the magnetic field strength, (right) The distribution of the H2O maser features indicating the 
velocity structure of the maser features. The open symbols are the red-shifted features and the solid symbols are the blue-shifted 
features. The size indicates the velocity difference with the stellar velocity (v ra d = 5.3 km s _1 ) 



clination and position angle of the magnetic field and the field 
strength at the surface of VX Sgr. The stellar radiu s was fixed at 
~ 16 A U, i n agreement wi t h the observations bv lDanchi et alJ 
ll 19941) and lCireenhill et alJ ( 1995). In the fits we also included 
the maser features for which we only determined an absolute 
upper limit. The errors on the magnetic field strength were 
taken from our analysis, while we included errors of 1 km/s 
on the velocity of the maser features to take into account turbu- 
lent velocities and small deviations from the spherical outflow 
model. Figs.[5]show the structure of the magnetic field as well 
as the radial velocity structure observed for the H2O maser fea- 
tures. The best fitted model for a dipole magnetic field is over- 
plotted. While the magnetic field strengths can be fit reasonably 
well, the model has a harder time simultaneously fitting the ve- 
locity structure. Our fits show the stellar position to be offset 
from the center of the H2O maser emission by * 70 x 70 mas to 
the NE. The magnetic axis of the dipole field is pointed toward 
us at an inclination angle i = 40 + 5°, and a position angle of 
= 220 + 10°. The indicated errors are the formal l<x errors 
from the fit. However, as both the magnetic field and the out- 
flow velocity structure are likely to have small scale structure, 
the exact errors are strongly model dependent and are hard to 
quantify. 

The fitted values are remarkably consistent with the val- 
ues found for the magnetic field determined from OH masers, 
as well as with the orientation angl e determined from the 
H2O maser distribution in M03 and lMarvell (l996). We find 
that the magnetic field strength at the surface of VX Sgr corre- 
sponds toBs; 2.0 + 0.5 kG. For the r~ 3 magnetic field strength 



dependence of the dipole field, this is consistent with the obser- 
vations of t he magnetic field on the OH masers as w ell as the 
SiO masers (Barvainis. Mci ntosh. & Predmorelll987i) . In addi- 
tion to a dipole magnetic field, we also considered a solar-type 
magnetic field with a magnetic field strength dependence of r~ 2 
on the distance to the star. This resulted in a significantly lower 
likelihood. 

Our ability to perform such a fit of VX Sgr indicates that 
an ordered large scale magnetic field exists. Previously, an 
alternative model explained the magnetism in terms of local 
fields, frozen in high density pockets. Although such small 
scale structure could still be important, we conclude that it is 
likely that large scale magnetic fields permeate the CSEs. 

6. Discussion 

6.1 . Circular polarization of H 2 masers 

The observations reveal circular polarization percentages up to 
~ 20%. Using the non-LTE radiative transfer method this in- 
dicates magnetic field strengths along the maser line of sight 
up to ~ 4 G, while the fields determined using the LTE ap- 
proximation are generally ~ 40% higher. As shown in V02, 
the non-LTE circular polarization spectra are typically not sym- 
metric, while the LTE method produces spectra with an point 
symmetric S-curve shape. However, due to the data process- 
ing and the necessary removal of a scaled down replica of 
the total power, it is impossible to directly observe the non- 
symmetric non-LTE spectra. Still, we confirm that the non- 
LTE spectra produce better fits to the observed V-spectra, as 



8 



W.H.T. Vlemmings et al.: Magnetic Fields around Late-type stars 




MilliARC SEC 
Center al RA 16 25 47.4772000 DEC 18 53 32.946000 
Peak flux = 9.8038E+00 JY/BEAM 
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Fig. 1. Total intensity image of the H2O maser features around 
U Her. The inlay shows the total intensity spectrum of the 
U Her H2O masers, where the solid vertical line indicates the 
stellar velocity and the dashed vertical line indicates the veloc- 
ity of the maser fea ture that was found to be coincid ent with 
stellar position in Iviemmings. van Langevelde. & Diamond! 
ll2002l) . 



can be seen in several of the spectra presented in Figs. |2] [5] 
and|4] Most of the observed V-spectra also show narrowing 
that cannot be reproduced in the LTE analysi s, as was first 
found in VOL Similar to earlier results (e.g. Spe ncer et all 
Il979l Iviemmings & van Langeveldel 120051) . we find that the 
H2O masers are mostly unsaturated after an analysis of the line 
widths and shapes of the total intensity spectra. 

We confirm in these observations again the absence of lin- 
ear polarization to a limit of * 2%, which provides compelling 
evidence against the non -Zeeman interpretation, considered by 
IWiebe & Watsonl Jl998l) for the SiO masers, as the cause for 
the observed high circular polarization of the H2O masers. The 
fact that no significant linear polarization is observed is an ad- 
ditional indication that t he circums tellar H2O masers are unsat- 
urated JNedoluha & Watsonl fl99ll) . 

6.2. The strength and shape of the magnetic field 

The observed magnetic field strengths on the H2O masers are 
consistent with the results for the other maser species assum- 
ing a B oc r" dependence of the field strength on the distance 
to the star. As the magnetic field strength determined for the 
SiO, H2O and OH maser features depends on the angle be- 
tween the magnetic field and the maser propagation axis, it is 
difficult to exactly determine the value of a. In Fig.|6]we show 




H„0 masers 



OH masers 



14 15 
log(R) [cm] 
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Fig. 6. Magnetic field strength, B, as function of distance, R, 
from the center of the star. Dashed-dotted boxes are the esti- 
mates for the magnetic fields in the OH and SiO maser regions 
of Mira stars, solid boxes are those for the supergiants. The 
dashed-dotted line indicates a solar-type magnetic field and the 
dotted line indicates the dipole field. The solid symbols are 
the previously observed stars (U Her: triangles; S Per: squares; 
VY CMa: crosses and NML Cyg: hexagonals). The open sym- 
bols indicate the stars of the sample presented here (U Her: tri- 
angle; VX Sgr: hexagonal and U Ori: circle). Note that the open 
triangle only indicates the upper limit determined for the most 
recent observation of U Her. Also note that the magnetic field 
strength observed on the OH masers of U Ori is larger than the 
typically observed fields for Mira variables. The dashed lines 
represent estimates of the stellar radius. 



an updated version of Fig. 15 from V02 including the results 
from the sources observed for this paper. The open symbols 
indicate the stars in the sample presented here. The points for 
the H2O masers are drawn at the outer radius of the maser re- 
gion, which for VX Sgr was determ ined in M03. For U Ori 
we used the value of « 30 AU from Bains et alJ J2003). The 
arrows indicate a typical width of the H2O maser region. The 
SiO maser location of x 50 AU for VX Sgr was taken from 
iDoeleman. Lonsdale. & Gr eenhillld l998l) and its OH maser ex- 
tent 1400 AU) from lSzvmczak Cohen. & Richardl d200ll) 
The O H mase rs around U O r i were found to exist at « 65 AU 
dChaoman. Cohen. & Saikial ll 99 ll) . For the other sources the 
references are listed in the relevant section of V02. We indi- 
cate a solar-type magnetic field (a = -2) and a dipole field 
(a - -3). The observations of VX Sgr indicate (as shown in 
§ |3 that a dipole field best fits the shape of the field in the 
H2O maser region. Extrapolating the observed magnetic field 
strengths to the stellar surface, we find that Mira variable stars 
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have surface field strengths up to several times 10 2 G, while 
supergiant stars have fields of the order of 10 3 G. 



6.3. The origin of the magnetic field 



The origin of the strong, large scale magnetic fields around 
evolved stars remains a topic of debate. The generation of an 
axisymmetric magnetic field requires a magnetic dynamo in 
the interior of the star. One of the main arguments against a 
dynamo generated magnetic field is the slow rotation of the 
late-type star based on conservation of ang ular momentum 
dSoker & Harpazlll992l) . However, models bv lBlackman et all 
(2001) consider the strong differential rotation created when 
the degenerated stellar core contracts while the envelope ex- 
pands. In this case the dynamo is the result of the interaction 
between the differential rotation (<y) and turbulence in the con- 
vection zone around the degenerated core (a). This aw-dynamo 
action can produce a strong magnetic field similar to the one 
observed using the circumstellar masers. It has been argued that 
such strong magnetic fields should produce st rong X-ray emis- 
sion. Although Mira is a weak X-ray source dSoker & Kastnerl 
120031) . observatio ns of TX Cam and T C as did not reveal any 
X-ray emission dKastner & S oker. 2004). However, as the op- 
tical depth for X-ray emission around evolved stars is expected 
to be high, it has been shown that this non-detection cannot 
rule out the strong magn etic fields produced by the <m>-dynamo 
(Blackma n et allEpQll) . The above models include the interac- 
tion with a degenerated core, while the core of supergiants such 
as VX Sgr is supposedly not degenerate. However, a similar dy- 
namo action for supergiants, driven by the differential rotation 
between the contracting non-degenerate core and the expand- 
ing outer envelope has been sho wn to also be a ble to generate 
strong magnetic fields llUchida & BaPDulll982l) . 

In contrast to the <m>-dynamo, several papers have exam- 
ined a o^w-dynamo, under the assumption that the rotation of 
the e nvelope only marginally contribute s to the dynamo ac- 
tion ( So ker & ZoabiL 120021 ISoken. 120021) . These models pro- 
duce magnetic fields that are several orders of magnitude lower 
than the models using the <m>-dynamo. The observed magnetic 
field strength in the maser regions of the CSE are then argued to 
be due to localized, magnetized wind clumps that are produced 
in magnetic spots on the surface of the star that can have local 
field strengths of up to ~ 10 G. Using another model, dynamo 
action of giant-cell convection at the surface of a late-type su- 
pergiant star has recently been shown to be able to produce 
fields u p to ~ 500 G in localized spots on the stellar surface 
(Dorch, 2004). Models that produce local strong field however, 
do not readily explain the observations of large scale magnetic 
field structure such as seen around VX Sgr. Alternatively, in the 
case of the <? 2 <jj-dynamo a strong magnetic field can be gener- 
ated when the star is spun up by a close binary. The sources 
in our sample do not however, show any indication of binarity, 
although this cannot be ruled out. 



6.4. Shaping of the circumstellar outflow 



The magnetic field strength of several Gauss measured in the 
H2O maser region implies that the magnetic pressure domi- 
nates the thermal pressure of the circumstellar gas throughout 
a large part of the CSE of both regular AGB stars and super- 
giants. The ratio between the thermal and magnetic pressure 
is given by /? = (SnnukT)/ B 2 ), with k the Boltzman constant. 
Assuming a gas density of nu - I0 g cirT 3 and a temperature of 
T « 10 3 K at the inner edge of the H2O maser region, a mag- 
netic field between B = 0.5 and 1 G gives B ~ 0.012 - 0.003, 
indicating that the magnetic pressure dominates the thermal 
pressure by factors of « 80 or more. In the lower density non- 
masing regions this factor will be even higher unless the mag- 
netic field is frozen into the high density clum ps. Still, in that 
case, when B oc n k with 1/3 < k < 1/2 ( Mousch oviasl ll 987l) . B 
will be of the same order of magnitude. 

The effects of magnetic fields on the stellar outflow and the 
shaping of the distinctl y aspherical PNe have been discussed 
in several papers (e.g . IPasco hi Il987t Ichevalier & Luol Il994l 
lGarcfa-Seguralll997l) . PNe are thought to be formed due to the 
interaction of the slow AGB wind (v ~ 10 km s _1 ) with a subse- 
quent fast wind (v ~ 1000 km s _1 ) generat ed when the ce ntral 
star e volves into a white dwarf dKwok. Purton. & Fitzgeraldl 
Il978h . Strong, large scale, magnetic fields around AGB stars 
coul d directly affect the fas t wind and as a result shape the 
PN dChevalier & Luol fl994l Idarcfa-Segura. L6pez.& Franco! 
2003). Only recently, magnetic field strengths of 1-3 kG have 
been detected at the surface of the central star of several 
PNe djordan. Werner & O'Toolet 120051) . Another mechanism 
in which mag netism can play a role in shaping PNe is described 
in Matt et all d2000l) . In the models presented in that paper, a 
large scale dipole magnetic field is responsible for creating an 
equatorial density enhancement in the initial, slow, AGB wind. 
For the magnetic field to influence the stellar outflow, the out- 
flow must be sufficiently ionized, which is likely the case close 
to the surface of the star. A dipole magnetic field will then only 
need to have sufficient magnetic pressure (B < 2.0) to produce 
an equatorial disk. The interaction of the later fast wind with 
this disk has been shown to be ab le to easily create aspheri- 
cal, cylindrical symmetric PNe (e.g . llckelll988HSoker & Livid. 
1989; Mellema, Eulderink, & Icke, 1991). Furthermore, under 
the influence of the large scale magnetic field, the circumstellar 
disk could become warped dLaii 120031) . and the i nteraction of 
the fast wind with such a disk has been shown bv llckd (2003) 
to be able to explain the multipolar shape observed for several 
PNe, such as NGC 7027 dCox et allEooli) . Although we have 
not directly observed the morphology of the magnetic field for 
the AGB stars in our sample, the magnetic field strengths mea- 
sured on U Ori, are fully consistent with a dipole magnetic field 
such as found around the supergiant VX Sgr. As seen in Fig|3] 
we observe indeed that the H2O masers around VX Sgr occur 
in a oblate spheroid, which could be an indication of the equa- 
torial density enhancement expected in a CSE that has been 
shaped by a dipole magnetic field. 
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7. Conclusions 

We have measured the magnetic field around the late-type stars 
U Her, U Ori and VX Sgr using observations of the circular po- 
larization of their H2O masers. Although we only find an upper 
limit in the case of U Her, we find strong magnetic fields of 
~ 3.5 and 4 G for U Ori and VX Sgr respectively. The rich 
structure of the H2O masers around the supergiant VX Sgr en- 
abled us to determine the shape of the magnetic field around 
this supergiant star. The observations are best represented by 
a dipole magnetic field at angles that are remarkably consis- 
tent with those of the dipole field used to explain previous 
OH maser polarization observations at much larger distances 
from the central star. This confirms the presence of an ordered 
magnetic field close to the star. Additionally, the ellipsoidal 
structure of the H2O masers around VX Sgr is aligned with 
the equatorial plane of the dipole field, which could indicate 
the equatorial density enha ncem e nt cau sed by the magnetic 
field as described by M att et al. (2000). The magnetic field 
strengths determined in the H2O maser regions around U Ori 
and U Her, lower mass evolved stars that are the progenitors 
of PNe, are also consistent with a dipole field, such as found 
around VX Sgr. As a result, also for those stars, the magnetic 
fields can cause aspherical density structures that result in non 
spherically symmetric PNe. 
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